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ABSTRACT 

We present the current estimate of instrumental and systematic effect uncertainties for the Planck-Low Frequency Instrument relevant to the first 
release of the Planck cosmological results. We discuss an overview of the main effects and of the tools and methods applied to assess the residuals 



in 



on maps and power spectra. We also present an overall budget of systematic effect uncertainties, which are dominated by far sidelobe straylight 
pick-up and imperfect calibration. They are at least two orders of magnitude less than the cosmic microwave background (CMB) temperature 
CO power spectrum. A residual spurious signal above the noise level is present, especially at 30 GHz, in the multipole range i < 20, and is likely to be 
caused by residual Galactic straylight contamination. Current analysis aims to further reduce the level of spurious signals in the data and to better 
CO model systematic effects, in particular straylight and calibration uncertainties. 
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1. Introduction 

This paper, one of a set associated with the 2013 release of data 
from the Planck^ mission (Planck Collaboration I 2013), de- 
scribes the Planck-LFl instrument systematic effects and their 
related uncertainties in CMB temperature maps and power spec- 
tra. 

The LFI implements a pseudo-correlation differential design 
similar to WMAP (Jarosik et al. 2003; Jarosik et al. 2003) to 
suppress l/f amplifier gain and noise fluctuations (Seiffert et al. 
2002; Mennella et al. 2003; BersanelH et al. 2010) as well as 
correlated effects from thermal and electrical variations affect- 
ing both the sky signal and reference loads. The reference signal 
is provided by stable 4.5 K blackbodies thermally and mechani- 
cally connected to the external structure of the High Frequency 
Instrument (HFI) 4 K box (Valenziano et al. 2009; Lamarre et al. 
2010). The offset between the sky and reference signals, of the 
order of 1 - 2 K, is balanced in software during data processing 
on the ground (Mennella et al. 2003; Zacchei et al. 2011). The 
differenced timestreams are characterised by 1 // noise knee fre- 
quencies in the range 10- lOOmHz (Mennella et al. 2010, 2011), 
leaving residual correlated low-frequency fluctuations in gain 
and signal that are removed during calibration and map-making. 

The LFI is also an excellent polarimeter, with very low sys- 
tematic effects. Depolarization by the optics and by imperfec- 
tions in the orthomode transducers, separating the orthogonal 
linear polarizations, has been accurately measured on the ground 
and is almost negligible (Leahy et al. 2010). 

Asymmetrical bandpass response in the two radiometers is 
the main source of / (Q,U) leakage in the foreground- 
dominated sky regions, especially at low frequencies. Although 
accurate knowledge of the bandpass response allows us, in prin- 
ciple, correction for this effect during data analysis, the ground 
bandpass measurements were not accurate enough to maintain 
this residual below 1% (Zonca et al. 2009). For this reason 
the spurious polarization from bandpass mismatch was esti- 
mated and removed using flight data, as described in Planck 
Collaboration II (2013). 

Optical effects arise mainly from Galactic and CMB dipole 
pick-up caused by primary and secondary mirror spillovers 
(Tauber et al. 2010; Sandri et al. 2010). This is relevant espe- 
cially for polarization measurements at 30 GHz, where Galactic 
emissions are stronger. 

In this paper we provide a comprehensive overview of the 
instrument systematic effects and the uncertainties they cause 

^ Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
entific consortia funded by ESA member states (in particular the lead 
countries France and Italy), with contributions from NASA (USA) and 
telescope reflectors provided by a collaboration between ESA and a sci- 
entific consortium led and funded by Denmark. 



on CMB temperature maps and power spectra (see Sect. 2). In 
Sect. 3 we outline and discuss the known instrumental effects, 
separating them into two broad categories: (i) effects that do not 
depend on the sky signal and impact the radiometric measure- 
ments as an additive spurious fluctuation or a gain variation, and 
(ii) effects that depend on the sky signal, i.e., on its amplitude 
and/or on the scanned sky region. Some of these effects are re- 
moved in the data processing pipeline according to algorithms 
described in Planck Collaboration II (2013). The assessment of 
the residual uncertainty, discussed in Sect. 4, was performed ac- 
cording to two different strategies. Null tests were the primary 
tool to check for systematic effect residuals exceeding the white 
noise level. We also assessed their impact on radiometric time- 
lines, even if below the white noise limit, by exploiting in-flight 
housekeeping and scientific data. This latter approach, that com- 
bines flight data with ground-measured instrument properties, 
has never been applied in previous CMB experiments and proved 
a powerful tool to check for systematic effects. 

Some of the effects discussed in this paper are also rel- 
evant for calibration, and are discussed in detail in Planck 
Collaboration V (2013). In this case we provide here only a brief 
discussion of the most relevant points and results, deferring to 
the dedicated paper any further details. 

Throughout this paper we follow the naming conven- 
tion described in Appendix A of Mennella et al. (2010) and 
also available on-line on the Explanatory Supplement (Planck 
CoUaboration ES 2013). 



2. Summary of uncertainties due to systematic 
effects 

In this section we provide a top-level overview of the uncertain- 
ties due to systematic effects in the Planck-LFl CMB temper- 
ature maps and power spectra. Table 1 provides a list of these 
effects, with short descriptions of their cause, strategies for their 
removal and references to sections and/or papers where more in- 
formation can be found. 

The impact of l/f noise has been assessed using "half-ring" 
noise maps (see Sect. 4.1.1) normalized to the white noise esti- 
mate at each pixel obtained from the white noise covariance ma- 
trix, so that a perfectly white noise map would be Gaussian and 
isotropic with unit variance. Deviations from unity trace the con- 
tribution of residual 1 // noise in the final maps, which ranges 
from 0.06% at 70 GHz to 2% at 30 GHz, as detailed in Sect. 12.2 
of Planck CoUaboration II (2013). 

Pixel uncertainties due to other systematic effects have been 
calculated on simulated maps degraded to N^idQ = 128 at 30 and 
44 GHz and A^side = 256 at 70 GHz in order to approximate the 
optical beam size. 

In Table 2 we list the r.m.s. and the difference between the 
99% and the 1% quantiles in the pixel value distributions. For 
simplicity we refer to this difference as the peak-to-peak (p-p) 
difference, although it neglects outliers but effectively approxi- 
mates the peak-to-peak variation of the effect on the map. 

The last row in Table 2 reports the r.m.s. value of survey 
difference maps (first minus second survey), after a correction 
which takes into account the white noise (see Sect. 4.1.2 for de- 
tails about this correction). These values thus represent the level 
of r.m.s. uncertainty per pixel due to known and unknown sys- 
tematic effects. 

Angular power spectra have been obtained from full reso- 
lution (A^side = 1024) systematic effect maps at each frequency 
using the HEALPix Anafast routine (Gorski et al. 2005). We 
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Table 1. List of known instrumental systematic effects in Planck-hFl 



Effect 



Source 



Control/Removal 



Reference 



White noise correlation .... 

1// noise 

Bias fluctuations 

Thermal fluctuations 

1 Hz spikes 

Sidelobe pickup 

Bandpass asymmetries .... 

Analogue-to-digital converter 
non linearity 

Imperfect photometric .... 
calibration 

Pointing 



Effects independent of sky signal 

Phase switch imbalance Diode weighing 

RF amplifiers Pseudo-correlation and destriping 

RF amplifiers, back-end electronics Pseudo-correlation and destriping 

4 K, 20 K and 300 K thermal stages Calibration, destriping 

Back-end electronics Template fitting and removal .... 



Effects dependent on the sky signal 



Main and sub-reflector spillovers . 

Difl'erential orthomode transducer 
and receiver bandpass response 

Back-end analogue-to-digital . . . 
converter 

Sidelobe pickup, radiometer noise 
temperature changes and other 
non-idealities 

Uncertainties in pointing reconstru- 
ction, thermal changes afl'ecting 
focal plane geometry 



Model sidelobes removed from timelines 
(not implemented in this release) 

Spurious polarisation removal 

Template fitting and removal 



Calibration using the 4 K reference 
load voltage output 

Negligible impact on temperature . 
anisotropy measurements 



3.1.1 
3.1.1 
3.1.3 

3.1.2, 4.2.1 
3.1.4, 4.2.3 

3.2.1,4.3.1 

Planck Collaboration II (2013) 
3.2.2, 4.3.2 

Planck Collaboration II (2013) 
3.2.3,4.3.3, 

Planck Collaboration V (2013) 
3.3,4.4 



Table 2. Summary of systematic eflTects uncertainties on maps^ 
in //KcMB- 



30 GHz 44 GHz 70 GHz 





p-p 


rms 


p-p 


rms 


P-P 


rms 


Bias fluctuations 


0.08 


0.01 


0.10 


0.02 


0.23 


0.06 


Thermal fluctuations . . 


. 0.61 


0.11 


0.40 


0.08 


1.17 


0.20 


1-Hz spikes 


0.87 


0.17 


0.14 


0.03 


0.60 


0.12 


Sidelobes pickup . . . . 


18.95 


4.53 


1.92 


0.57 


6.39 


1.91 


ADC non-linearity . . . 


. 3.87 


1.01 


0.89 


0.19 


0.92 


0.19 


Gain residuals 


4.33 


1.16 


4.74 


0.97 


6.51 


1.10 


Total^ 


21.02 


4.83 


5.61 


1.13 


7.87 


2.00 


Survey consistency^ . . 




8.50 




7.10 




8.20 



^ Calculated on a pixel size approximately equal to the average beam 
FWHM. 



^ The total has been computed on maps resulting from the sum of indi- 
vidual systematic eff'ect maps. 

^ Measured r.m.s. of the half-diff'erence between two complete sky sur- 
veys minus the r.m.s. of the noise for pixels, see Sect. 4.1.2. 



have then evaluated the propagation of the various effects in the 
final CMB map by assuming a simple internal linear combina- 
tion component separation, as explained in Sect. 4.5. In Fig. 1 
we show how the power spectra of the various eflTects compare 
with the Planck temperature spectrum, with the noise level com- 
ing from the half-ring diflTerence maps (see Sect. 4. 1 . 1) and with 
the residual map obtained from a diflference map between survey 



1 and survey 2^ (see Sect. 4.1.2). The large plot in the top panel 
shows the power spectra obtained from frequency-independent 
maps resulting from the weighted- average of frequency maps us- 
ing the weights specified in Sect. 4.5. Spectra in the three small 
plots in the lower panel, instead, show contributions of system- 
atic eflTects from individual frequency maps. 

Our assessment shows that the global impact of systematic 
eflfect uncertainties is at least two order of magnitudes less than 
the CMB power spectrum, and demonstrates the robustness of 
Planck-LFl temperature anisotropy measurements. Comparison 
between the total simulated systematic eflTects and residual sig- 
nal obtained by diflTerencing survey 1 and survey 2 maps high- 
lights an excess signal in the multipole range £ < 20 that is not 
completely accounted for in our simulations. This excess comes 
mainly from the 30 GHz channel and is likely to be caused by 
Galactic emissions picked up by beam sidelobes. Understanding 
this excess and further reducing the level of residual systematic 
uncertainties is the primary goal of our current analysis to obtain 
polarization measurements with a level of purity comparable to 
what has been achieved for temperature anisotropics. 



3. Overview of LFI systematic effects 

Known systematic effects in the Planck-LFl data can be divided 
into two broad categories: effects independent of the sky signal, 
which can be considered as additive or multiplicative spurious 
contributions to the measured timelines, and effects which are 
dependent on the sky and that cannot be considered indepen- 
dently of the observational strategy. 



^ Time periods relative to individual surveys are defined in Table 1 1 
of Planck Collaboration II (2013) 
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00 




CMB -TT 


— (E) Bias fluct. 




— (A) Gain uncertainty 


— (F) Temp, fluct. (4 K) 


— (SD) Survey difference 


— (B) Far sidelobes 


— (G) 1-Hz spikes 




— (C) ADC 


— (H) Temp, fluct. (300 K) 


— (HR) Half-ring noise 


— (D) Temp, fluct. (20 K) 


Total 





Fig. 1. Angular power spectra of the various systematic effects compared to the Planck temperature anisotropy spectrum. The 
black dashed curve, representing the total contribution, has been derived from a map where all the systematic effects have been 
summed. Top panel: power spectra obtained from frequency independent maps resulting from the weighted- average of individual 
systematic effect frequency maps. Bottom panel, contributions of systematic effects from individual frequency maps. The CMB 
curve corresponds to the Planck best- fit model presented in Planck Collaboration XV (2013). 
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Fig. 2. Amplitude spectral density estimates on 5 -day time pe- 
riods (coloured lines) compared to the nominal mission noise 
model for a representative 70 GHz radiometer (LFI23M). 
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Fig. 3. Amplitude spectral density estimates on 5 -day time pe- 
riods (coloured lines) compared with the nominal mission noise 
model (black line) for one 44 GHz radiometer (LFI24S). 



3.1. Effects independent of sky signal 
3.1 .1 . Noise correlations and 1// noise 

Each Planck-hFl receiver is a pseudo-correlation system view- 
ing a scalar feed directed through the telescope at the sky, to- 
gether with a reference cold load thermally stable near 4 K. Non- 
white noise from the cold front-end amplifiers is reduced via 
the correlation, while fluctuations in the later stages of the re- 
ceiver are minimized by modulating a phase switch in the corre- 
lation section at 8192 Hz. The LFI receiver design, construction, 
ground performance and initial flight performance have been ex- 
tensively documented (Bersanelli et al. 2010; Mennella et al. 
2010, 2011). 

The noise properties of the receivers play an important role 
in downstream data analysis. In particular, we need good es- 
timates of the white noise level, long term stability (l//-type 
noise) and any correlated noise components. 

The receiver architecture is symmetric, with two comple- 
mentary detector diodes as output for each receiver channel. As 
described in Seifl'ert et al. (2002) and Mennella et al. (201 1) im- 
perfect matching of components generates isolation between the 
complementary diodes of a receiver between -10 and -15 dB. 
This imperfect isolation leads to a small anti-correlated compo- 
nent in the white noise that is cancelled by a weighted average of 
the time ordered data from the two diodes of each receiver as the 
first step of analysis. This avoids the complication of tracking 
the anti correlated white noise throughout the analysis. 

We treat the combined diode data as the raw data, and cal- 
ibration, noise estimation, map-making etc. are performed on 
these combined data. The weights were determined from some 
initial estimates of the calibrated noise for each detector, and are 
kept fixed for the entire mission. 

Noise parameters were reported in Mennella et al. (201 1). A 
longer data set, some thermal instabilities in the instrument (par- 
ticularly during survey 3), and refinements of the data analysis 
(map making and noise co variance matrix) all require a more de- 
tailed look at the long term evolution of the noise characteristics 
of the receivers. 



The noise power spectral density P(f) of the receivers is gen- 
erally well described by 



P(f) = o-^ 



(1) 



where cr characterizes the white noise component, the knee fre- 
quency, /k, denotes the frequency where white noise and l/f 
contribute equally in power to the total noise, and a character- 
izes the slope of the power spectrum for frequencies f < fk- 

We estimate the signal- subtracted noise power spectrum of 
each receiver on 5-day time periods. Except for specific, mostly 
well understood events, shorter time scale noise estimation does 
not produce any evident trends. For nearly all the radiometers 
our noise model is a very good approximation of the power spec- 
trum. We plot a representative comparison in Fig. 2. A few chan- 
nels show features not well captured by this simple model; the 
worst is displayed in Fig. 3. 

Over the course of the nominal mission, the noise is well fit 
by the model, with the exception of the early parts of sky sur- 
vey 3. During this time, thermal instabilities brought on by the 
switch-over from the nominal to the redundant sorption cooler 
cause poor fits and some changes in the parameters. In Figs. 4 
through 6 we show the behaviour of the three noise parameters 
in Eq. (1) estimated on 5-day sections of data over the nominal 
time period. White noise and knee frequency are stable, while 
the slope starts increasing in absolute value after day 300, as a 
result of larger temperature fluctuations in the 20 K focal plane. 
The jump in slope after day 500 is correlated with the sorption 
cooler switch-over (see Sect. 3.1.2 for further details). 

3.1.2. Thermal effects 

The LFI is susceptible to temperature fluctuations in the 300 K 
back-end modules, in the 4 K reference loads and in the 20 K 
focal plane. Figure 7 provides an overview of the main temper- 
atures during the period between day 91 (the start of nominal 
operations) and 563 after launch. 

The two topmost plots show the reference load tempera- 
tures at the level of the 70 GHz and 30 -44 GHz channels, re- 
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Fig. 4. Fitted white noise parameters over the nominal survey 
for representative radiometers at 30, 44 and 70 GHz. Values are 
estimated on 5 -day sections of data. 
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Fig. 5. Fitted knee frequencies over the nominal survey for rep- 
resentative radiometers at 30, 44 and 70 GHz. Values are esti- 
mated on 5 -day sections of data. 



spectively. The temperature of the 70 GHz reference loads is ac- 
tively controlled by a proportional-integral-derivative (PID) sys- 
tem and is very stable (ST^jris ~ 0.13mK, see the zoomed plot 
in the inset). Reference loads of the 30 and 44 GHz channels, 
instead, do not benefit from active thermal control. Their tem- 
perature is consequently more unstable and susceptible to major 
system-level events like, for example, the switchover to the re- 
dundant sorption cooler. 

The third plot from the top of Fig. 7 shows the 20 K LFI fo- 
cal plane temperature measured by a sensor placed on the feed 
horn flange of the LFI 2 8 receiver. The temperature during the 
first sky survey was very stable, with a ^rnns ^ 1 mK. Towards 
the end of the first year of operations the sorption cooler perfor- 
mance started to degrade and its stability was maintained with a 
series of controlled temperature changes. The switchover to the 
redundant cooler was performed on August, 1 1* 2010, leaving a 
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Fig. 6. Fitted power-law slopes for low frequency noise. Here 
we note significant instability during the early part of the third 
survey. This is due to substantially greater thermal instability of 
the 20 K stage just before and after switchover between the two 
sorption coolers. 



clear signature on all the main LFI temperatures. After this op- 
eration the level of temperature fluctuations in the focal plane 
increased unexpectedly, and this was later understood to be the 
efl'ect of liquid hydrogen that was still present in the cold-end 
of the nominal cooler, because the degraded compressor system 
was not able to absorb all the hydrogen that was present in the 
cooler line. Although this efl'ect was later mitigated by a series of 
dedicated operations, most of the third sky survey sufl'ered from 
a higher- than-nominal level of temperature variation. 

The last plot shows the temperature of the 300 K electron- 
ics box, measured by one of its temperature sensors. During the 
first sky survey the back-end temperature suff'ered from a daily 
fluctuation caused by the satellite transponder that was switched 
on daily during contact with the ground station. After day 258 
the system was left continuously on, and the modulation dis- 
appeared. This operation caused an increase of the absolute tem- 
perature level. The second temperature change occurred in corre- 
spondence to the sorption cooler switch-over operation. The plot 
also shows a yearly temperature modulation due to the satellite 
rotation around the Sun and a temperature spike at day 191 after 
launch. This was caused by an operational anomaly that led the 
satellite to fail to repoint for an entire day with a corresponding 
temperature increase of the warm units. 

More details about the thermal stability performance of 
Planck can be found in Planck Collaboration II (201 1), while the 
susceptibility of the LFI to temperature variations is discussed in 
Terenzi et al. (2009). 

3.1.3. Bias fluctuations 

The signal detected by the radiometers can vary because of fluc- 
tuations in the front-end and back-end amplifier bias voltages. In 
the LFI these fluctuations occurred according to two time scales: 

- slow electric drifts, due to thermal changes in the power sup- 
ply, in the RF amplifiers, and in the detector diodes; 

- fast and sudden electric instabilities, arising in the warm 
electronics or from electromagnetic interference eff'ects, and 
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Fig. 7. Main temperatures in Planck-LFl. From top to bottom: 
70 GHz reference loads, 30 and 44 GHz reference loads, 20 K 
focal plane (sensor placed on feed horn flange of LFI 2 8) and 
300 K back-end (sensor placed on the back-end electronics box). 
A brief description of the main operational events afl'ecting the 
thermal behaviour is provided in each panel. 



aff'ecting both the cold amplifiers and the warm detector 
diodes. 

The eff'ect of slow drifts is suppressed by the pseudo- 
correlation architecture of the diff'erential radiometers. Fast elec- 
tric changes produce quasi-random fluctuations and abrupt steep 
drops or jumps in the signal. If jumps are caused by instabilities 
in the front-end bias voltage then the efl'ect involves the output 
voltage of both diodes in the radiometer. When the jumps occur 
in the back-end detector diodes (so-called "popcorn noise") they 
impact only the output voltage of the corresponding diode and 
afl'ect sky and reference load samples. In both cases the difl'er- 
enced signal is largely immune from these efl'ects. 

3.1.4. 1-Hz spikes 

This efl'ect is caused by pickup from the housekeeping elec- 
tronics clock that occurs after the detector diodes and before 
the analogue-to-digital converter (ADC) (Meinhold et al. 2009; 
Mennella et al. 2010, 2011). This spurious signal is detected in 
the radiometer time-domain outputs as a 1 s square wave with a 
rising edge near 0.5 s and a falling edge near 0.75 s in on-board 
time. In the frequency domain it appears at multiples of 1 Hz. 

Frequency spikes are present at some level in the output from 
all detectors, but afl'ect the 44 GHz data most strongly because 



of the low voltage output and high post-detection gain values 
in that channel. For this reason spikes are removed from the 
44 GHz time-ordered data via template fitting, as described in 
Planck CoUaboration II (2013). 



3.2. Effects dependent on sky signal 
3.2.1 . Sidelobe pick-up 

Stray light contamination arises from the spurious signal pickup 
from the telescope far sidelobes. Main sources of stray light con- 
tamination are the Galaxy, especially at 30 GHz, and the cos- 
mological dipole, mainly detected in the directions of the main 
and sub-reflector spillover, as sketched in Fig. 8. In principle 
we should also include the straylight contribution from the or- 
bital dipole, but its efl'ect is a factor ten lower than the cosmic 
dipole, so that it can safely be neglected in this framework (but 
it has been considered in the calibration pipeline). Further de- 
tails about the Planck optical system are reported in Tauber et al. 
(2010), while the LFI beam properties are provided in Sandri 
et al. (2010). 



Main reflector 
spillover 



Sub-reflector 
spillover 




Fig. 8. Main and sub-reflector spillover, and main beam direc- 
tions in the Planck telescope. 

Straylight impacts the measured signal in two ways: (i) 
through direct contamination and coupling with the main beam 
sky signal, and (ii) in the photometric calibration of the radiome- 
ter detected signal. In this paper we concentrate on the direct 
detection, while the impact on calibration and the adopted miti- 
gation strategies are described in Planck Collaboration V (2013). 

Because of the beam orientation, the straylight fingerprint is 
diff'erent in odd surveys compared to even surveys. The Galaxy, 
for example, is detected by the sub-reflector spillover in the odd 
surveys and by the main-reflector spillover in the even surveys. 
Because the sub-reflector spillover points approximately in the 
main beam direction, the Galaxy straylight pattern is close to 
the Galactic plane. The main-reflector spillover, instead, points 
at about 85° from the main beam so that the Galaxy is re-imaged 
onto a ring (see figures in Sect. 4.3.1). 

3.2.2. ADC non linearity 

The ADC linearity requires that the voltage step sizes between 
successive binary outputs are constant over the entire input dy- 
namic range. If these steps are not constant (see the sketch in 
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Fig. 9) we have a non linearity in the ADC response that leads to 
calibration errors. A brief description of the mathematical model 
of this effect is provided in Appendix A. 




Voltage output 



Fig. 9. Schematic of the ADC non-linearity effect. For a small 
range of voltages the ADC response changes slope. 

The typical fingerprint of ADC non-linearity is a variation of 
the detector voltage output white noise not paired by a detectable 
variation in the voltage level. This effect was observed in the 
LFI radiometer data for the first time in flight, where drops of 
a few percent were observed in the voltage white noise but not 
in the output level over periods of few weeks. Fig. 10 shows 
this effect as a plot of relative white noise variation versus the 
detector output voltage for one the the most affected radiometer 
channels (the 44 GHz detector LFI25M-01). 

The grey points represent an average over each pointing 
(about 40 min) while the solid line has been obtained by further 
binning the data in 200 bins over the plotted range in order to 
reduce the scatter and show more detail. The figure shows that 
the typical amplitude of the region where the non-linearity oc- 
curs is of the order of 1 mV, corresponding to about three bits in 
the ADC. The ADC effect is strongest (3 to 6%) in the 44 GHz 
channels, because of their lower detector voltages. 
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Fig. 10. Plots of the percentage variation in the single detector 
white noise estimates with detector voltage. 



The ADC non-linearity effect has been characterised from 
flight data and removed from the data streams according to the 
procedure described in Planck Collaboration II (2013). In Fig. 1 1 
we show the same data as in Fig. 10 after the correction has been 
applied. The figure clearly shows that the anomalous white noise 
dips disappear after correction. 

3.2.3. Imperfect photometric calibration 

An important set of systematic effects are those related to the 
photometric calibration of the radiometers. Such effects are dis- 
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Fig. 11. Same as in Fig. 10 after correction of the ADC non- 
linearity effect. 



cussed in Planck Collaboration V (2013); here we will only pro- 
vide the most important information to put the results of that 
paper within the context of this work. 

There are three different kinds of systematic effects that can 
affect the calibration. 

1. Incorrect assumptions regarding the calibration signal. In 
the case of LFI, the signal used for the calibration is the dipo- 
lar field caused by the motion of the Solar System with re- 
spect to the CMB rest frame and by the motion of the space- 
craft around the Sun. We model the former using the val- 
ues quoted by Hinshaw et al. (2009) and the latter using the 
spacecraft's attitude information. Any error in the numbers 
would directly lead to an error in the calibration of Planck- 
LFI data. 

2. Incorrect treatment of the calibration signal. To actually use 
any previous knowledge of the CMB dipole, we need to con- 
volve the signal with the beam response of the LFI radiome- 
ters. Any error in this step would produce a systematic ef- 
fect in the map, not only because of the wrong shape ex- 
pected for the calibration signal, but also because of the re- 
moval of the (wrong) dipole from the calibrated maps done 
by the Planck-LFl pipeline (Planck Collaboration II 2013). 
Possible types of errors include: systematic errors in WMAP 
dipole estimate, wrong convolution of the expected dipole 
with the radiometer beams, incorrect masking of the Galaxy 
when fitting the observed signal with the dipole, etc. 

3. Incorrect reconstruction of gain fluctuations . Some of the 
algorithms we used in calibrating LFI data for this release 
use the radiometer equation and the recorded variations of 
the radiometers total-power output to track gain changes. In 
principle, any deviation in the behaviour of the radiometer 
from the ideal case (e.g., ADC non-linearities) can induce 
systematic effects in the gain curves. 

3.3. Pointing effects 

Pointing uncertainties are translated into uncertainties in pixel 
temperature measurements. If pointing uncertainties are not con- 
stant in time then the statistics of the sky anisotropy measure- 
ments is not preserved, with a consequent impact on power spec- 
trum and cosmological parameters. For Planck-LFl, pointing un- 
certainties arise from two main effects: 

1. Satellite pointing determination. The Planck Attitude 
Control Movement System guarantees a pointing accuracy 
of about (Planck Collaboration I 2013; Planck Science 
Office 2010), which is well within scientific requirements. 
However, small non-idealities in the system and errors in 
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the attitude reconstruction (caused, for example, by thermo- 
elastic effects) can affect the data. 
2. Uncertainties in the focal plane geometry reconstruction. 
The measurement of the Planck-LFl focal plane geometry is 
based on the determination of the beam pointing with respect 
to the nominal line of sight exploiting Jupiter observations. 
The peak of each beam has been determined by fitting data 
with a bivariate Gaussian function which may not be repre- 
sentative of the real beam centre. 



white noise covariance maps (Planck Collaboration ES 2013) 
and checking the standard deviation of the resulting maps. We 
found this r.m.s. value to be very close to unity: 1.021 1, 1.0089, 
and 1.0007 for 30, 44, and 70 GHz, respectively. The deviation 
from unity is consistent with the different level of l/f noise in 
the three frequency channels (see Tables 1 and 10 of Planck 
Collaboration II 2013). 



4. Assessing residual systematic effect 
uncertainties in maps and power spectra 

In this section we discuss the assessment of the impact of resid- 
ual systematic effects on maps and power spectra. This assess- 
ment has been performed according to two strategies: by "null 
maps" obtained by differencing maps with the same sky signal, 
in order to highlight residuals, and by simulating the systematic 
effects in the timelines, exploiting a combination of flight data 
and measured instrumental properties. 

4.1. Null tests 

We define a "null test" as any difference between two indepen- 
dent data sets which are anticipated to give nearly the same sig- 
nal, on the assumption of perfect calibration, pointing recon- 
struction, and systematic effects removal. Null tests are a pow- 
erful means to assess the validity and self-consistency of Planck 
data on various timescales and across different dimensions (de- 
tector, frequency, time), and to highlight systematic effects above 
the white noise level. 

The design of Planck and its observing strategy provides a 
wide range of opportunities for null tests, with sensitivity to dif- 
ferent systematic effects and implications for the scientific out- 
puts. Although we refer to these tests as null tests, the results 
are generally not featureless, and are required to match our in- 
strument model, in terms of far sidelobe contribution, main lobe 
ellipticity and orientation, and especially our noise model, cap- 
tured in the noise covariance matrix produced with each map. 

In the Planck-LFl collaboration each internal data release is 
accompanied by a comprehensive set of null tests as a check of 
our processes and ongoing improvement in terms of systematic 
errors. In this section we report the results from the main tests 
supporting the systematic effect analysis for the first Planck pub- 
lic data release. 

Here we provide some sample difference maps. Unless other- 
wise noted, these maps are masked to remove point sources and 
to include only pixels measured in both maps. Difference maps 
are divided by 2 to be statistically consistent with average maps, 
and are smoothed to 2° FWHM to enhance large scale features. 

4.1 .1 . Half-ring difference null tests 

Half-ring difference null tests, constructed by taking a weighted 
difference between the first and second halves of each pointing 
period, are useful to assess the data noise properties and system- 
atic effects on time scales smaller than about 20 min. Weights 
are calculated as explained in Sect. 9.2 of Planck Collaboration 
II (2013). 

In Fig. 12 we show the half-ring difference maps for the 
three LFI frequencies. Visual inspection of these maps shows 
no apparent correlated structures. A quantitative test was per- 
formed by dividing them pixel-by-pixel by the square root of the 
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Fig. 12. Half-ring difference maps: 30 GHz (top), 44 GHz (mid- 
dle), and 70 GHz (bottom). 



4.1 .2. Survey difference null tests 

Differences between single survey maps are useful to check for 
residual systematic effects at large angular scales. 

Odd minus even surveys difference maps highlight effects 
arising from beam ellipticity and far sidelobes, as discussed in 
Sect. 3.2.1. In Fig. 13 we show the survey 1 minus survey 2 dif- 
ference maps for the three LFI frequencies. The cosmic and or- 
bital dipole signals are removed during calibration (as discussed 
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in Planck Collaboration V 2013), so the difference in orbital 
dipole signal between survey 1 and survey 2 is not visible. Far 
sidelobe pickup of the galactic plane in survey 2 is visible as a 
large blue ring which matches the simulation results shown in 
Figs. 24 through 26. 
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Fig. 14. Survey 1-2 difference map at 30 GHz minus the survey 
1-2 difference map obtained from 30 GHz simulated sidelobes 
(Fig. 24). The blue ring disappears. 



known and unknown systematic effects after correcting for the 
white noise contribution (Xnoise- The latter has been estimated us- 
ing the high-^ tail of the survey half-difference power spectra, by 
means of the following formula: 

noise - ^ , K^) 

where A^pix is the number of pixels in the map and (Q) is the av- 
erage level of the power spectrum in the white-noise dominated 
range 2000 < £ < 3000. This value of (Tnoise quantifies the r.m.s. 
due to statistical noise in the survey difference maps, and it is 
equal to^ 14.0, 16.2, and 28.1yuK/pixel at the frequencies of 30, 
44, and 70 GHz, respectively. We then use the following formula 
to estimate the r.m.s. due to systematic effects only (Table 2), 




-10 mKcmb 10 



Fig. 13. Survey 1 minus survey 2 difference maps: 30 GHz (top), 
44 GHz (middle), and 70 GHz (bottom). Such odd-even differ- 
ences highlight beam ellipticity and sidelobe pickup effects. 



Figure 14 shows the result of a test performed on the 30 GHz 
survey 1 minus survey 2 null test (first map of Fig. 13). The 
test consisted in removing from this map a difference map be- 
tween survey 1 and survey 2 simulated sidelobes (see Fig. 24 
in Sect. 4.3.1). The result (reported also in Planck Collaboration 
IV 2013) showed that the blue ring structure disappears, which 
confirms both the nature of this spurious feature in the map and 
that our sidelobe model is correct with a good approximation. 

From the survey difference maps at the three LFI frequency 
bands we have estimated the r.m.s. uncertainty per pixel due to 



^syst ~ ~ ^noise' 

where cr sdis is the r.m.s. of the survey difference maps. The equa- 
tion holds exactly only if all the signals are gaussian, but numer- 
ical simulations have shown that Eq. (3) is true within 10-15 % 
for the kind of effects considered in this work. 

In Fig. 15 se show survey 1 minus survey 3 null test maps. 
Survey 1 and survey 3 cover the sky in nearly identical orien- 
tations, and would be consistent with noise if calibration and 
other systematics were perfectly controlled. The visible large 
scale features appear to be due to the relatively poor thermal 
stability of the instrument for the first portion of survey 3. 

We can also use angular power spectra to estimate the im- 
pact of residuals seen in our difference maps. The spectrum of 
the half-ring difference maps provide an accurate estimate of 
the noise contribution. A range of systematic effects we might 
be concerned about, such as imperfectly tracked gain changes, 
systematic thermal variations and others, might be expected to 
change significantly over the course of 1 year of operations. 
These spectra show that, to very good approximation, the signal 
in the survey difference maps (basically data taken six months 
apart) is fully explained by the instrumental noise. The fact that 
the half-ring difference spectra and survey difference spectra are 
nearly identical gives us confidence that neither the systematic 
effects we know about, or others we may not have anticipated, 
are contaminating the data in any significant way. 

^ The values reported here use A^pix = 12 A^^^.^^, with A^side equal to 128 
(30, 44 GHz) or 256 (70 GHz). 
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Fig. 16. Angular power spectra for 30 GHz null tests. Pseudo- 
spectra are calculated on 80% of the sky, with the Galaxy and 
point sources masked. All spectra are corrected for sky fraction 
but not for beam smearing effects. A best fit cosmological model 
spectrum is displayed for scale. For ^ > 30 the spectrum of the 
survey 1 minus survey 2 map diff'erence lies directly on top of 
the half-ring diff'erence spectrum calculated for the same time 
period, the (absolute value) diff'erence of these spectra is also 
shown. For this 30 GHz spectrum, the lower £s are contaminated 
by the sidelobe leakage at the 20 juK^ level. 
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Fig. 15. Survey 1 minus survey 3 diff'erence maps: 30 GHz (top), 
44 GHz (middle), and 70 GHz (bottom). 



4.2. Assessment of timeline-additive systematic effects 
4.2.1 . Thermal effects 

Method. Thermal systematic effects maps have been gener- 
ated using a simulation strategy that combines in-flight tempera- 
ture sensor measurements (Mennella et al. 2010), thermal mod- 
elling of the propagation of temperature fluctuations (Tomasi 
et al. 2010) and radiometric transfer functions measured during 
ground tests (Terenzi et al. 2009). Here we sketch the procedure 
used to combine these data into systematic effect maps. 

For each temperature effect and for each receiver detector 
diode we have chosen the most representative sensor, generally 
the closest to the receiver. Housekeeping data have been low- 
pass filtered to remove high frequency sensor noise, Fourier- 
filtered to obtain the estimated temperature fluctuation at the re- 
ceiver location, and then multiplied by the radiometric transfer 




Fig. 17. Angular power spectra for 44 GHz null tests. For this 
frequency, the low £ spectrum of the survey 1 minus survey 2 
map difference is closer to that of the half-ring difference spec- 
trum. Letters in the plot follow the same convention of the legend 
in Fig. 16. 



function to obtain the simulated antenna temperature fluctuation 
on the undifferenced sky and reference load channels. 

For each pointing period the average measured sky and refer- 
ence load voltages have been added to the two antenna temper- 
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Fig. 18. Angular power spectra for 70 GHz null tests. At 70 GHz, 
the survey 1 minus survey 2 map difference and half-ring differ- 
ence spectra are essentially the same everywhere. Letters in the 
plot follow the same convention of the legend in Fig. 16. 



ature fluctuation data streams. After a weighted average of the 
two detector data values of each radiometer, we took the sky- 
load difference using the gain modulation factor, r, and multi- 
plied the resulting stream by the photometric constant, G. The 
weights, r, and G, were the same as used in the nominal pipeline 
to produce sky maps for that radiometer. 

After oversampling to the receiver sampling frequency using 
linear interpolation, we used these data to build maps using the 
in-flight pointings and the same map making procedure used to 
produce the final scientific products. 



□ 300 K back end 

□ 4 K reference loads 




Fig. 19. Peak-to-peak thermal effects on maps. Notice the loga- 
rithmic scale on the ordinate axis. 
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Results. In Fig. 19 we show the peak-to-peak amplitude of the 
various effects on final maps. 

Back-end temperature fluctuations have a sub-yuK effect on 
maps. This low level can be understood if we consider that these 
fluctuations impact sky and reference load signals symmetrically 
and are effectively suppressed in the differential measurement. 
Furthermore the residual present in the data acts as a purely mul- 
tiplicative effect, so that it is essentially calibrated out by apply- 
ing our gain model (Planck Collaboration V 2013). 

Temperature variation at the level of the 4 K reference loads 
couple with the radiometric output as an asymmetrical additive 
spurious signal. In this case the relative calibration model pro- 
vides no benefit, so that at 30 and 44 GHz this effect leaves a 
residual of about 1 yuK peak-to-peak. The effect at 70 GHz, in- 
stead, is largely suppressed by the active thermal control system 
present on the HFI focal plane close to the reference loads of this 
frequency channel. Front-end 20 K temperature variations cou- 
ple with the radiometric measurements through both gain and 
noise temperature fluctuations. For this reason only a part of the 
effect can be calibrated out. Moreover, the asymmetry of the re- 
ceiver chain before the orthomode transducer is such that the 
effect suppression provided by the sky-load differencing is not 
optimal. The residual effect is similar for the three frequency 
channels and is of the order of 1 jjK peak-to-peak. 

Maps of the combined thermal effects at the three LFI fre- 
quency channels are shown in Fig. 20. 
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Fig. 20. Maps of combined thermal effects at 30 GHz (top), 
44 GHz (middle), and 70 GHz (bottom). 
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4.2.2. Bias fluctuations 



Method. The effect of bias fluctuations in the front-end ampli- 
fiers has been computed on maps and power spectra using the 
measured drain currents and a linear transfer function that links 
the drain currents of the two amplifiers to the radiometric out- 
put in antenna temperature. Since we are interested in assessing 
purely electrical instabilities, we have corrected the drain current 
housekeeping data to remove variations induced by temperature 
changes in the 20 K and 300 K temperature stages, i.e.. 



CTn(0 = /drain(0 " C^20K^r20K(0 " C^SOOK^^SOOkCO, (4) 

where ^r20K(0 and ^r3ooK(0 are temperature variations on the 
20 K and 300 K temperature units, respectively, and 0^20 k and 
0^300 K are the corresponding drain current thermal susceptibil- 
ity coefficients, which were calculated using an iterative linear 
fitting process. First we calculated the coefficients of the sus- 
ceptibility to back-end temperature fluctuations, exploiting the 
temperature change induced by the change in the transponder 
state, which occurred at day 258 (see Fig. 7), and then we deter- 
mined the coefficients of the susceptibility to front-end tempera- 
ture fluctuations using data from a temperature susceptibility test 
run at the end of the in-flight calibration phase. The iterative pro- 
cess was closed by re-calculating back-end thermal coefficients 
after correcting drain currents for front-end temperature fluctua- 
tions. 

After thermal correction we correlated drain current chang es 
with antenna temperature variations in sky and reference load 
samples. We recall here that the LFI receiver architecture im- 
plies that the signal characteristics at each detector depend on 
both the radiometer front-end amplifiers (Bersanelli et al. 2010). 
So, for each detector diode we have first calculated the weight, 
w, providing the maximum correlation between the output volt- 
age and the linear combination of the drain currents of the 
two radiometer amplifiers, then a linear fit between this com- 
bination and the voltage output provided the required transfer 
function. Mathematically the relationship between the corrected 
drain current fluctuations (^^^^ain^ voltage output varia- 

tions (^VskyCref)) rcadS! 
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^Vsky(ref)(0 = ^sky(ref) [w 61^^ ,0) ^ (k - w) 61^^ ,0)] , (5) 

where ^ is a constant, w is the weight and asky(ref) is the slope of 
the linear fit between the weighted combination of the two drain 
currents and the sky (reference load) voltage outputs. The time- 
ordered data obtained by Eq. (5) have then been projected onto 
the sky using flight pointings. 

Results. Maps of the residual effect arising from bias fluctua- 
tions are shown in Fig. 21. This residual effect is less than 1 jjK 
peak-to-peak at all frequencies and presents little structure apart 
from a stripe in the 70 GHz map, that is caused by a jump in the 
bias voltage that occurred at day 258, following the change in 
transponder state. This jump affected more the 70 GHz radiome- 
ters so that a small signature was left in the maps. 

4.2.3. 1-Hz spikes 

Method. Time ordered data containing the spike signal have 
been generated using templates obtained from flight radiometric 
data. The details of this method have been described in Sect. 7.1 
of Mennella et al. (201 1) and will not be repeated here. 



-0.1 /iKcMB 0.1 

Fig. 21. Maps of the effect from drain current fluctuations at 
30 GHz (top), 44 GHz (middle) and 70 GHz (bottom). 

Results. In Fig. 22 we show maps of the spike systematic effect 
at the three LFI frequencies. Because, as described in Planck 
Collaboration II (2013), spikes are removed from 44 GHz but 
not from 30 and 70 GHz time-ordered data, the corresponding 
maps represent the residual effect after removal at 44 GHz and 
the spike effect with no removal applied at 30 and 70 GHz. In all 
the three channels the r.m.s. effect is at the sub-yuK level. 

4.3. Assessment of effects dependent on the sky 
4.3.1. Sidelobes 

Method. The external stray light contamination has been evalu- 
ated with simulations in which the sky model includes the dif- 
fuse Galactic emission and the dipole, which are the most impor- 
tant sources of external straylight contamination. At 30 GHz, the 
straylight contamination has been assessed including the beam 
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32, and 33 GHz). The optical model used in GRASP simulations 
is the "RFTM" reported in Planck Collaboration IV (2013). 

For each frequency interval an observation of the sky model 
was simulated for all 30 GHz detectors using the beam sidelobes 
and the real sky pointings, neglecting beam smearing effects and 
weighting the data stream with the above mentioned weight fac- 
tors. Finally the Madam map-making code was run to generate 
maps from simulated data streams. 
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Fig. 23. Bandpass response of the two radiometers of the LFI 2 7 
receiver. The figure shows the seven frequency intervals and the 
corresponding frequencies at which sidelobes have been simu- 
lated. For each interval the weight is the integral of the bandpass 
response curve. 
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Fig. 22. Maps of 1 Hz spikes at 30 GHz (top), 44 GHz (middle) 
and 70 GHz (bottom). The map at 44 GHz represents the residual 
after the spike signal has been removed from the time-ordered 
data, while maps at 30 and 70 GHz represent the spike signal 
with no removal applied. 

frequency dependence and the receiver in-band response. The 
bandpass response (Zonca et al. 2009) of the radiometers was di- 
vided into discrete frequency intervals. For each frequency inter- 
val a weight factor was calculated as the integral of the bandpass 
response over the interval itself. In Fig. 23 we show, as an exam- 
ple, the bandpass response of the LFI 2 7 receiver (main and side 
arms) and the seven frequency intervals considered in the simu- 
lations. The weights correspond to the integral of the bandpass 
response curve over the frequency interval. 

In parallel, far sidelobes were computed using the GRASP 
MrGTD"^ software (www . ticra . com) at the frequencies indicated 
on the top of each slice reported in Fig. 23 (27, 28, 29, 30, 31, 

^ Multi-reflector geometrical theory of diffraction 



Results. In Fig. 24 through 26 we show the simulated sidelobe 
fingerprint on the sky, after the destriping process, for the odd 
surveys (left side) and even surveys (right side), respectively. 
The figures show that the straylight from the cosmological dipole 
is similar in the two surveys, while the Galaxy straylight is larger 
in the second survey, as expected. The ring-shaped fingerprint in 
the second survey has also been revealed at the expected level 
in the real data by taking the diff'erence of even minus odd sur- 
vey maps (see Figs. 13 and 14 in Sect. 4.1), thus confirming the 
reliability of our simulations. 




-7.4 //KcMB 7.4 -7.4 /xKcmb 7.4 



Fig. 24. 30 GHz channel: sidelobe fingerprint on the sky due to 
Galactic foregrounds (top row) and cosmological dipole (bottom 
row). Survey 1 (left side) and survey 2 (right side). 

Results show that the most sensitive channel to straylight is 
30 GHz, followed by 70 GHz and 44 GHz, which is the most in- 
sensitive channel to straylight eff'ects. This is consistent with the 
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Fig. 25. 44 GHz channel: sidelobes fingerprint on the sky due to 
Galactic foregrounds (top row) and cosmological dipole (bottom 
row). Survey 1 (left side) and survey 2 (right side). 
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Fig. 26. 70 GHz channel: sidelobes fingerprint on the sky due to 
Galactic foregrounds (top row) and cosmological dipole (bottom 
row). Survey 1 (left side) and survey 2 (right side). 

telescope optical performance at the various frequencies. In fact, 
the primary mirror is strongly under-illuminated by the 44 GHz 
horns, resulting in a low stray light sensitivity at the expense of 
a larger main beam, especially for the LFI25 and LFI26 horns. 
The 30 and 70 GHz horns are characterised by similar illumina- 
tion properties, so that their stray light susceptibilities are com- 
parable, with a slightly better performance of the 70 GHz horns 
with respect to the 30 GHz ones. If we also take into account the 
larger sensitivity of the 30 GHz channel to the Galactic signal it 
is apparent that this channel is, overall, the most susceptible to 
stray light contamination. 

Finally, to quantify the straylight eff'ect on maps and power 
spectra, we have generated a global map per frequency, includ- 
ing the dipole and Galactic straylight signals for both surveys 
(see Fig. 27). 

4.3.2. ADC non-linearity 

Method. The levels of the residuals due to ADC correction pro- 
cess were estimated by applying the correction algorithm to sim- 
ulated data containing a known ADC eff'ect and making diff'er- 
ence maps with those produced from data with no ADC eff'ect. 
Ring-based sky and reference load time-ordered data for all in- 
dividual detectors were constructed. The noise component was 
estimated using average 1// noise parameters over the nominal 
mission for each detector. Galactic and CMB signals were rep- 



systematic uncertainties 
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Fig. 27. StrayHght contamination maps at 30 GHz (top), 44 GHz 
(middle) and 70 GHz (bottom). 



resented by using the observed maps, which were rolled out in 
the time-ordered data using real pointings and successively un- 
calibrated using the inverse of the appropriate gain table. The 
same was done with a map of the WMAP dipole, while the or- 
bital dipole was calculated from the pointing information and 
JPL ephemeris for the satellite velocity. Finally spline fits to 
the observed sky and reference voltage levels per pointing pe- 
riod were used together with estimated receiver temperatures, 
^CMB = 2.725 K and T^cf = 4.5 K, as a model for the gain evolu- 
tion. 

An ADC eff'ect was induced by applying the inverse of the 
spline correction used in the real data. The same algorithm as 
used with the real data was applied to the simulated data iter- 
atively five times to ensure convergence. Intensity maps were 
constructed just using a simple binning of the data into an 
A^side = 1024 map at each iteration and also to data with no 
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ADC effect added. Some of these maps showed a residual dipole 
caused by small changes in the overall slope of the temperature- 
voltage response curve due ADC correction. Since the calibra- 
tion pipeline will determine this response and not give rise to a 
residual dipole, we also remove it here via a correlation fit with 
the input dipole map. The ADC eff'ect maps are then the diff'er- 
ence between the fifth iteration map and the no- ADC map. Maps 
for each frequency band are then produced by averaging all maps 
for that frequency, taking into account the detector weighting. 

Since not all the 70 GHz channels could be corrected due 
to some channels having "popcorn" noise, another method was 
used to estimate the likely level of ADC error, using the white 
noise level on the diff'erence data. This is immune to the "pop- 
corn" noise, but cannot be used to correct the ADC eff'ect, since 
it is not known whether the eff'ect is on the sky or reference volt- 
ages. In this case we only estimate the ADC eff'ect and do not 
apply any correction. 

Results. Maps of the ADC eff'ect at the three LFI frequency 
channels are shown in Fig. 28. The main eff'ect of the ADC 
residuals is a small (< 0.1% of the dipole signal) ring-based 
gain error which appears in the maps as stripes in the scan di- 
rection. The residuals are generally larger where the strongest 
signals are: the dipole and Galactic plane. The contribution from 
the Galactic plane becomes weaker at higher frequencies as ex- 
pected. Broad stripes seen in the 30 GHz map are due to residual 
deviations from linearity on voltages ranges larger than the ADC 
peaks. These also occur at the other frequencies, but as the num- 
ber of channels increases this eff'ect averages out, leaving more 
uniform noise-limited, low-level residuals at 70 GHz. While the 
44 GHz channels have the strongest ADC eff'ect due to lower de- 
tector voltages, they are also the best characterized, leading to a 
well-determined correction placing it between 30 and 70 GHz in 
terms of the amplitude of residuals. 

4.3.3. Imperfect photometric calibration 

Method. We have developed an analytical model of the im- 
pact of the uncertainty in the dipole calibration algorithm due 
to the radiometer white noise and the loss of integration time 
due to Galactic masking. We have run this model to estimate 
how this effect propagates through the calibration and mapmak- 
ing pipeline. Such simulations scan a sky map (the input map) 
of pure astrophysical signal (without dipole) to produce a time- 
ordered data stream, which is then uncalibrated using gains in- 
ferred from the total-power output of the radiometers. These 
time-ordered data are then used as input in a simplified version 
of the LFI pipeline to produce a new calibrated map (the output 
map). The difference between the input and output maps should 
be mainly due to dipole leakage, since the gains used in the de- 
calibration phase differ from those calculated by the pipeline. 
Refer to Sect. 5 in (Planck Collaboration V 2013) for more in- 
formation. 



Results. Fig. 29 shows the difference between the input and 
output maps. The shape of the features in these maps closely 
follows the scanning circles drawn by the pointing direction of 
the telescope towards the sky. (This is expected, since the cal- 
ibration is performed on the time-ordered data.) The estimated 
impact of such systematic effects on the Planck-LFl maps is of 
the order of a few yuK per pixel. 




-0.5 /iKcMB 0.5 

Fig. 28. Maps of the ADC non-Hnearity effect at 30 GHz (top), 
44 GHz (middle), and 70 GHz (bottom). 

4.4. Pointing uncertainties 

Method. To estimate the uncertainty introduced by the 
Gaussian fit in main beam measurements we performed tests 
using the radio-frequency model of the flight telescope (Planck 
Collaboration IV 2013) and we compared the centre calculated 
by the fit with the beam maximum, which is uniquely determined 
in optical simulations. Typical differences between the centres 
are V for afl the 70 GHz beams, 4'' for LFI24, 18'' for LFI25 
and LFI26 (44 GHz horns), and 6'' for the beams at 30 GHz. All 
these estimates are less than the statistical uncertainty in the de- 
termination of the beam centre, which ranges from 4'' at 70 GHz 
to 10'' at 30 GHz. 

The focal plane geometry was reconstructed using four 
Jupiter transits labelled as Jl, J2, J3 and J4 (Planck Collaboration 
II 2013; Planck CoUaboration IV 2013). When we compared the 
focal plane geometry obtained from the combination of Jl and 
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1. generate time ordered data by observing a CMB-only sky 
with flight detector pointing derived by applying the two fo- 
cal plane database solution; 

2. reconstruct the CMB map from the time-ordered data gener- 
ated in step 1 applying the two focal plane database solution 
in map reconstruction; 

3. repeat step 2 using the single focal plane database solution; 

4. take the difl'erence of the power spectra obtained from the 
two generated maps. 



-2.3 
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-2.3 
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X, 




-3.6 



2.9 



Fig. 29. Maps of the eff'ect of calibration uncertainties at 30 GHz 
(top), 44 GHz (middle), and 70 GHz (bottom). 



Results. Figure 30 shows that the relative difl'erence of power 
spectra obtained from the two maps generated according the 
above procedure is of the order of 10"^, which is negligible. 
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Fig. 30. Relative difl'erence between the spectra of the maps sim- 
ulated with one and two instrument databases. The three curves 
represent power spectra relative to three difl'erent feed horn pairs 
in the 70 GHz frequency channel. A running average smoothing 
has been applied to reduce the scatter and enhance any larger- 
scale trends. The order of magnitude of relative uncertainties is 
AQ/Q < 10-4. 



J2 with the one obtained from the combination of J3 and J4 we 
found a difl'erence of about 15'' in pointing, mainly along the 
in-scan direction. On the other hand, the comparison of the fo- 
cal plane geometries determined from single Jupiter transits (Jl 
against J2 and J3 against J4) showed difl'erences within the ex- 
pected uncertainty. The 15'' discrepancy, likely to be correlated 
with changes in the thermal control set-point of the data pro- 
cessing unit in the instrument digital electronics, was compen- 
sated using two difl'erent instrument databases in the data analy- 
sis pipeline, one for the period ranging from day 91 to day 539 
after launch, and the other for the period between day 540 and 
day 563. Details of the focal plane reconstruction and related 
uncertainties can be found in Planck Collaboration II (2013). 

We checked the impact of this efl'ect using dedicated simula- 
tions according to the following procedure: 



4.5. Propagation of systematic uncertainties in component 
separation 

A further step in our assessment has been to evaluate the im- 
pact of the various systematic efl'ects on the CMB map indepen- 
dently from the frequency. In order to do this we have weight- 
summed the three maps for each efl'ect using weights obtained 
using a pixel-based ILC (internal linear combination) compo- 
nent separation (Leach et al. 2008; Planck Collaboration XII 
2013). ILC relies on variance minimisation exploiting a multi- 
frequency dataset including a flat spectrum CMB component in 
thermodynamic temperature plus foregrounds characterized by a 
non-thermal spectrum. With this approach the CMB temperature 
can be determined at each pixel, p, by a simple weighted sum of 
the frequency maps once the correct weights are found: 
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Tcmb(p) = ^WiTy.(p), 

N 

The propagation of systematic effects in component separa- 
tion can be assessed using Eq. (6), in which we use the system- 
atic effect maps at the various frequencies instead of the mea- 
sured frequency maps. Therefore, for each pixel, we have that 
the residual systematic effect after component separation will be 
given by 

where the set of weights w/ is the same as used for retrieving 
the CMB map. In our assessment we have considered only LFI 
systematic effects and so we assumed, in Eq. (7), rsyst,y, = for 
Vi > 70 GHz. For the simplicity of this estimation, the weights 
are uniformly distributed in pixel and harmonic domains. In the 
Planck component separation pipeline, the variance minimisa- 
tion is conducted in the needlet space, i.e., on sub-sets of the 
harmonic and pixel domains where foregrounds are relevant at 
various levels, resulting in a set of coefficients for each needlet 
domain (Planck Collaboration XII 2013). 

4.6. Gaussianity statistical tests 

We now assess non-Gaussianity induced by the systematic ef- 
fects on the LFI data maps. We present results on different non- 
Gaussianity tests carried out on maps containing the various sys- 
tematic effects considered through the paper. In particular, we 
studied the map obtained by summing all the systematic effect 
maps. 

For detailed information on the non-directional or targeted 
non-Gaussianity tests on the Planck data see, for example, 
Planck Collaboration XXIII (2013), Planck Collaboration XXIV 
(2013), and Planck Collaboration XIX (2013). The Hst of es- 
timators for the non-Gaussianity tests is based on a selection of 
well-known statistics that are able to detect non-Gaussian signals 
with a wide range of amplitude and morphological properties. In 
particular we consider the Minkowski functionals (Schmalzing 
& Gorski 1998), statistical quantities derived from the 1 -point 
PDF (variance, skewness, kurtosis and the Kolmogorov-Smirnov 
or KS distance) and the skewness and kurtosis of the spher- 
ical Mexican hat wavelet (SMHW, Martinez-Gonzalez et al. 
2002). The properties of these estimators are described in Planck 
Collaboration XXIII (2013) and references therein. 

We have compared values of our estimators applied to a set 
of realistic Gaussian CMB and noise realizations and another set 
based on the same CMB and noise simulations summed with the 
systematic effect maps. The CMB and noise maps were simu- 
lated following the Planck-LFl data processing pipeline (Planck 
Collaboration II 2013). Where the three estimators did not pro- 
vide significant deviations between the maps with and without 
the systematic effects, we have carried out an additional test by 
reseating the systematic effect maps with a constant factor in or- 
der to provide an estimate of the needed amplitude to detect sig- 
nificant deviations with respect to the CMB signal (larger than 
3(T or 99% confidence level). 



Table 3. Impact of systematic effects on mean 1 -point PDF esti- 
mators 



Frequency [GHz] 


30 


44 


70 


Standard deviation 


-0.07 


-0.01 


-0.01 


Skewness 


-0.02 


-0.01 


-0.01 


Kurtosis 


-0.02 


0.00 


-0.01 



Figure 31 shows the three Minkowski functionals for the 
three LFI frequency bands. In each panel we compare the ±lcr 
(68%) confidence band centred on the mean corresponding to 
the Gaussian CMB plus noise simulations, with the same sim- 
ulations with systematic effects added for HEALPix resolutions 
A^side = 512, 256, and 128. Our analysis shows that the estima- 
tors based on the Minkowski functionals are not affected by the 
presence of systematic effects in the maps. 

Table 3 contains the difference of the mean of the two distri- 
butions (maps with and without systematic effects) normalized 
by its dispersion corresponding to the four statistics based on the 
1 -point PDF, namely standard deviation, skewness, kurtosis and 
KS distance for the three LFI frequency bands. The maps are at 
HEALPix resolution of A^side = 1024. There are no significant 
deviations, the distributions corresponding to the two types of 
map are almost superimposed and therefore the estimators based 
on the 1 -point PDF are not affected by systematic effects. 

Table 4 contains the difference of the mean of the two distri- 
butions (maps with and without systematic effects) normalized 
by its dispersion and corresponding to the skewness and the kur- 
tosis of the SMHW for the LFI frequency bands. The maps have 
a HE ALP i x resolution of A^side = 1 024 . The list of angular scales 
selected for this analysis is the one used in Planck Collaboration 
XXIV (2013) for the /ni estimation, and it is based on a set of 16 
angular scales between 1.3' and 956.3' logarithmically spaced. 
There are no significant deviations and therefore also the estima- 
tors based on the SMHW are not affected by systematic effects. 

Finally we have characterized the levels of detectability of 
the non-Gaussian contamination of these systematic effect maps. 
Using the statistics based on the 1 -point PDF and the Minkowski 
functionals, we have analysed simulations with different levels 
of systematic effects: 

Ar(n) = ArcMB(n) + Arno,se(n) + fAT,y,,(n), (8) 

where / is a factor that measures the level of detectability. We 
consider that the non-Gaussian signal is detected for the value of 
/ for which any of the estimators is outside the 3cr confidence 
level of the values corresponding to maps without systematic 
effects. The results indicate that the minimum value of the factor 
/ beyond which we detect non-Gaussianity is: / ~ 8 at 30 GHz, 
/ ~ 12 at 44 GHz, and / ~ 7 at 70 GHz. In summary, systematic 
effects do not generate significant levels of non-gaussianity for 
the temperature maps at the three LFI frequencies. 

5. Conclusions 

In this paper we have analysed and quantified the uncertainties 
on Planck-LFl CMB temperature anisotropy measurements aris- 
ing from systematic effects. The analysis has been fulfilled using 
two complementary approaches. 

Firstly we adopted a top-down approach, in which spurious 
excess signals were highlighted by a series of dedicated null- 
tests in which maps containing the same sky signal were dif- 
ferenced to obtain maps containing noise and systematic effect 
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Fig. 31. The three Minkowski functionals computed for Gaussian CMB and noise simulations (black symbols) compared with the 
Minkowski functionals computed for the same simulations with systematic effects added (solid magenta line). From left to right: 
the area, contour length or perimeter, and the genus. From top to bottom: the three LFI frequencies, 30, 44, and 70 GHz. These are 
explicitely for: A^side = 512. 



residuals. Secondly we followed a bottom-up approach in which 
each effect has been simulated to produce timelines and maps 
which have been further propagated through ILC component 
separation to a frequency-independent map. 

Our analysis has shown that systematic effect uncertainties 
are two orders of magnitudes or less below the CMB temper- 
ature anisotropy power spectrum, and are dominated by stray - 
light pick-up from far sidelobes and imperfect photometric cal- 
ibration. In this current data release the sidelobe signal is not 
removed from the data, but the CMB dipole pickup by far side- 
lobes is accounted for during the calibration process using a 
monochromatic model. 

Statistical analyses performed on maps containing the sum of 
all the simulated systematic effects added to a simulated CMB 
map showed no detectable non Gaussianity levels unless their 
level was artificially increased by a factor ranging from 7 to 12, 
which confirms that instrumental effects do not significantly im- 
pact Gaussianity studies. 

Comparison of power spectra obtained via the top-down and 
bottom-up approaches show a signal excess in the survey dif- 
ference maps at 30 GHz that is not accounted for in the simu- 



lated maps, in the multipole range i < 20. This excess could 
be caused by yet un-modelled straylight pick-up affecting the 
measurements both directly and in the photometric calibration 
process. 

Currently, analysis focusses on understanding and further re- 
ducing the level of systematic uncertainties in view of the 2014 
data release which will include polarization data and results. 
Areas of activity include a more thorough in-band modelling of 
the sidelobe response at all frequencies aimed at removing the 
spurious signal from timelines and at a better correction during 
the calibration step. 
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Table 4. Impact of systematic effects on skewness and kurtosis 
at various angular scales. 



Aungular scale 

n 




Skewness 






Kurtosis 




jyj v_rnz 


H-H- vjnz 


1 u vjnz 


j\j \jr\L 


A A ni-fy 


/ u vjnz 


1.3 


-0.13 


-0.20 


-0.06 


-0.39 


-0.40 


-0.25 


2.1 


-0.09 


0.09 


0.58 


0.08 


-0.42 


-0.10 


3.4 


-0.29 


0.13 


0.74 


0.61 


-0.04 


-0.06 


5.4 


-2.25 


-0.08 


0.99 


0.60 


0.04 


-0.88 


8.7 


-7.52 


0.05 


4.39 


1.30 


0.96 


-6.62 


13.9 


-6.52 


0.35 


1.86 


-0.20 


0.75 


-2.38 


22.3 


-1.23 


-0.05 


0.07 


-0.79 


-0.17 


-0.32 


35.6 


-0.20 


0.11 


-0.12 


-0.44 


-0.04 


-0.06 


J /.U 


n 1 n 


U.ZJ 


rv o 1 
— U.Zi 


— U.ZJ 


rv rvrv 


rv 1 ^ 
— U.io 


91.2 


0.19 


0.26 


0.02 


-0.02 


0.25 


-0.08 


146.0 


-0.07 


0.09 


0.14 


-0.13 


0.07 


0.07 


233.5 


-0.06 


-0.09 


0.33 


-1.34 


0.05 


-0.33 


373.6 


0.45 


0.28 


0.20 


-1.56 


0.08 


-0.21 


597.7 


0.80 


0.36 


0.16 


-0.01 


0.27 


-0.20 


956.3 


0.08 


0.21 


0.25 


-0.05 


-0.28 


0.12 



Values represent the normalized difference (in %) of the mean for the 
skewness and the kurtosis for each scale of the SMHW, considering 
maps with and without systematic effects. 



ities in which they have been involved, can be found at http : //www. sciops . 
esa . int/index . php?pro j ect=planck&page=Planck_Collaboration. 
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Appendix A: Theory of the ADC non-linearity effect 

The ADC non-linearity arises when the measured detector volt- 
age differs from the true voltage in some repeatable but com- 
plicated manner, depending on the exact positions of the voltage 
thresholds of the chip. By finding a way to map the apparent volt- 
age, V\ to the true voltage, V, the ADC eff'ect can be corrected 
and this mapping is just the ADC response curve, RiV'), as seen 
through the LFI acquisition system. In a perfect radiometer this 
voltage is also the product of the system temperature, Tsys, and 
radiometer gain, G(0, 

V=YR{Y) = G{t)T,y,. (A.l) 

Probing the response function requires tracking how a small 
known input voltage variation, AV, appears as a measured AV 
at diff'erent working voltages, V. The way this works can be 
illustrated by diff'erentiating the above equation with respect to 
V\ to obtain 



AV = I -h R{Y) I ^Y = G(t) AT. 



(A.2) 



Equation (A.2) shows the relation between diff'erential input 
and output signals and how a localized change in gradient can 
dominate via the dR/dV term. It also shows that small intrinsic 
thermal noise fluctuations, AT, can be used as the correspond- 
ing test input temperature signal assuming it to be due to band- 
width limited noise power Ar = T^y^/ VAv t, where Ay and r are 
channel bandwidth and sample integration time, respectively. By 
combining the two previous equations, the difl'erential response 
can be expressed as 



dR(r) 
dV 



Wavtav' y/ 

In the case of no ADC efl'ects and voltage variations induced 
purely through gain fluctuations, we have V = VAvt AV, and 
the difl'erential response dR(V^)/dV^ remains zero for all V\ as 
expected. Non-linearities are signaled where the thermal white 
noise does not follow detector voltages, revealing variations in 
the response curve. Since the radiometer gains drift very slowly, 
many estimates of detector white noise by Fourier analysis from 
the one minute scan rings are available, and by binning and aver- 
aging good signal-to-noises of 100 are achievable. The above 
equation can be integrated numerically, making use of these 
binned values, if the factor R(V') is taken to be unity, which is 
a good approximation. A discrete set of corrected voltages Vk 
for each of the binned measured voltages V- can be found via a 
trapezoidal summation as follows: 



Vk 



sr 



AV 



1 



1 



(A.4) 
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Here is lowest voltage bin, SV the voltage bin width, and 

a = 1/ VAvt is fitted so that the top voltage bin, Vmax, is equal 
to y^^x to maintain the same overall linear response. These ta- 
bles of corrected voltages to measured voltages for each detector 
are stored in the DPC database as the ADC correction, and are 
implemented as spline fits to these control points to carry out the 
correction on the time-ordered data. 
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